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We study the mean-field theory, and the properties of fluctuations, in an out of equilibrium 
Bose-Fermi system, across the transition to a quantum condensed phase. The system is driven out 
of equilibrium by coupling to multiple baths, which are not in equilibrium with each other, and 
thus drive a flux of particles through the system. We derive the self-consistency condition for an 
uniform condensed steady state. This condition can be compared both to the laser rate equation 
and to the Gross-Pitaevskii equation of an equilibrium condensate. We study fluctuations about the 
steady state, and discuss how the multiple baths interact to set the system's distribution function. 
In the condensed system, there is a soft phase (Bogoliubov, Goldstone) mode, diffusive at small 
momenta due to the presence of pump and decay, and we discuss how one may determine the field- 
field correlation functions properly including such soft phase modes. In the infinite system, the 
correlation functions differ both from the laser and from an equilibrium condensate; we discuss how 
in a finite system, the laser limit may be recovered. 
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I. INTRODUCTION 

In the last decade there have been enormous advances 
in the experimental realisation and theoretical under- 
standing of the phenomenon of quantum condensation, 
i.e macroscopic occupation of a single quantum mode, 
in different physical conditions. The phenomena ranges 
from Bosc-Einstcin Condensation (BEC) of structureless 
bosons to the BCS-typc collective state of fcrmions and 
has been studied in several physical systems such as de- 
generate atomic gases and supcrconductorai. Further, 
recent experimental advances in manipulation of atomic 
Fermi gases have led to realisation of the BCS-BEC 
crossover regime^ii^ and low-dimensional atomic conden- 
sates have also been explored^i^ii. From the early days 
of experimental investigation of BEC there have been 
enormous efforts in order to realise quantum condensa- 
tion in the solid stated. For this, the currently promising 
candidates are excitons in coupled quantum wells^i^, 
microcavity polariton o^^d^'^'^ , quantum Hall bilayersii, 
and Josephson junction arrays in microwave cavities^^. 
Although all these systems potentially may condense at 
temperatures orders of magnitude higher than those for 
dilute atomic gases, it has proven to be much more diffi- 
cult to realise BEC in the solid state than in atomic traps. 
Recently a comprehensive set of experiments^ reports 
polariton condensation in CdTe based microcavities but 
still the level of control in the study of the condensed 
states in solid-state is far from the finesse achieved in 
atomic vapours. 

In these various candidates for condensation, one 
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should distinguish different classes of systems. Equilib- 
rium superconductors are special in that the decay of 
pairs is disallowed. In equilibrium particlc-hole conden- 
sates, such as quantum Hall bilayers or charge density 
waves, particlc-hole mixing (tunnelling in bilayers) leads 
to a gapped spectrum; however the gap may be very 
small. Non-equilibrium particle-hole condensates in the 
solid state are, to a much greater extent than atomic 
gases, subject to dephasing and decay. It is not usually 
possible to isolate the condensate from the environment: 
lattice phonons, impurities and imperfections of the crys- 
tal structure lead to dephasing, and due to poor trapping, 
particles escape, requiring external pumping to sustain a 
steady-state. The dephasing and decay processes are of- 
ten faster than thermalisation. putting the system out 
of thermal equilibrium. The decay, and consequent lack 
of equilibrium have for a long time presented the ma- 
jor experimental obstacle in the realisation of solid-state 
condensation in otherwise appropriate conditions. Even 
if one can accelerate thermalisatio n^^i^^ , comparing the 
decay rates to other energy scales, one may see that de- 
cay, and the consequent flux of particles through the sys- 
tem, remains a more important effect in solid state than 
in atomic gases. 

Thus, a significant presence of dissipation and decay 
also poses fundamental questions about the robustness of 
a condensate, for example: whether a steady-state con- 
densate is possible with incoherent pumping and decay, 
and if so, how does it differ from thermal equilibrium, and 
from a lascr^iSi. Quantum condensation in dissipativc sys- 
tems also provides a connection to other phenomena of 
collective behaviour in the presence of dissipation such 
as pattern formatio n^^i^" , particularly in Iasers21ii2^i2^, 
and also recently in a system related to that studied 
here, the coherently pumped polariton optical paramet- 
ric oscillato r^^i^^ . Other recent examples of phase tran- 
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sitions and coherence in driven systems include quan- 
tum criticality in magnetic systems in the presence of 
current o^^'^^ , and transport through a Kondo dot coupled 
to multiple dots^. The relation between lasing and BEC 
is particularly relevant for polariton BEC, where the ex- 
perimental distinction between the two is not straight- 
forward^Si. 

Microcavity polaritons in particular, being made from 
fermionic particles and photons, have several special fea- 
tures and so provide an excellent laboratory to study 
condensation in dissipative environment. Due to the 
large wave-length of their photonic component and non- 
linearities associated with underlying fermionic structure 
the physics exits the regime of weakly interacting bosons 
at even modest density^. Putting aside a few subtleties 
characteristic only for polaritons one can say that with 
increasing density the quantum condensation transition 
moves from BEC (fluctuation dominated) to something 
like the BCS (mean-field, interaction dominated) col- 
lective stated, analogous to the BCS-BEC crossover in 
atomic Fermi gases near Feshbach resonance^i^. This al- 
lows one to explore the influence of non-equilibrium and 
dissipation not only on the usual BEC but also on more 
exotic forms of quantum condensation. A further compli- 
cation is that polaritons in planar microcavities are two- 
dimensional (2D) particles and so in an infinite equilib- 
rium system, although there is a Berezhinskii-Kosterlitz- 
Thoulcss (BKT) transition to a supcrfluid phase, below 
the transition long-wavelength fluctuations destroy the 
off-diagonal long range order and result in algebraic de- 
cay of phase coherence. Dissipation changes the struc- 
ture of collective modes and influences the spatial and 
temporal coherence in 2D quasicondensates, changing the 
power-law controlling the decay of phase correlationsi^. 
Finally, microcavity polaritons can also be trapped cither 
in stress-induced harmonic potentials^'^'^ or in natu- 
ral traps provided by microcavity disorder which reduce 
the influence of long- wavelength fluctuations and may al- 
low the existence of a true condensate and phase coher- 
ence over the whole system sizei^. How this confinement, 
when combined with pumping and decay, modifies the 
properties of coherence in such systems is an interesting 
questioii^, which has not yet been fully addressed. 

The last issue is particularly relevant for the deeper un- 
derstanding of the differences and connections between a 
polariton condensate and the laser. Apart from the obvi- 
ous difference; the laser being a collective coherent state 
of massless non-interacting photons while the condensate 
consists of massive and interacting bosons (in polariton 
condensation both massive photons and strongly coupled 
excitons are coherent), there are more subtle differences 
connected with fluctuations and so expected differences 
in the decay of correlationsi^. Lasing is normally consid- 
ered in systems with a well-defined single or a few mode 
structure and so the phase fluctuations which control the 
laser lincwidth are those of a phase diffusion of a sin- 
gle mode^. In contrast, condensation is usually studied 
in systems where there is a continuum of single parti- 



cle modes, and thus collective excitations involve coher- 
ent interaction of these different modes which affects the 
decay of coherence and the line-shape of the emissioni^. 
While lasing in systems with transverse freedom has been 
investigated for its pattern forming properties'^, there re- 
main many open questions concerning the decay of cor- 
relations and the crossover from a small system with few 
spatial modes to the infinite and many-mode limit. 

Although semiconductor microcavities in strong cou- 
pling provide a natural system to explore such phenom- 
ena, all these issues are by no means restricted to polari- 
ton condensation. With recent advances in manipulating 
dilute atomic gases similar conditions can be engineered, 
an immediate example is that of an atom laser in which 
a continuous leakage of atoms from atomic BEC takes 
place. To our knowledge the description of the output 
from an atom laser has been to date largely analogous 
to that of the photon laser— and the influence of the 
continuum of modes connected with atomic BEC on the 
coherence properties of the atom laser has not been ad- 
dressed. 

In a previous paper— we addressed some of these is- 
sues. We used a model Bose-Fermi system coupled to 
independent baths, not in thermal or chemical equilib- 
ria with each other, providing incoherent pumping and 
decay. We show that steady-state spontaneous conden- 
sation can occur in such systems, and can be distinct 
from lasing: The condensate can exist at low densities, 
far from the inversion required for lasing. We also found 
that the collective modes are qualitatively altered by the 
presence of pumping and decay: The low energy phase 
mode (Goldstone, Bogoliubov mode) becomes diffusive at 
small momenta. By considering the effect of phase fluc- 
tuations, we described the decay of correlations, which 
at large times and distances differs both from that for a 
thermal equilibrium condensate and from a laser. 

In this manuscript, apart from providing technical de- 
tails of the method we address several new aspects of 
quantum condensation in dissipative environment. In 
particular we study the influence of the cxciton density 
of states, and the temperature of the pumping bath on 
the non-equilibrium phase diagram. We also analyse how 
the non-thermal occupation of photon states is controlled 
by competition between the pumping and decay baths, 
and how this occupation deviates from that in thermal 
equilibrium. We do not a priori assume that the system 
is close to equilibrium, and so the system's distribution 
function may be of any form. Finally we provide a full 
account of how to determine field-field correlation func- 
tions in the condensed state, where phase fluctuations 
may be large, and so expansion to second order is in- 
sufficient. These field-field correlation functions describe 
the decay of correlations at large times and distances, 
and their Fourier transform gives the line-shape of a non- 
equilibrium condensate. This is an important extension 
to the non-equilibrium path integral techniques which to 
our knowledge has not been done before. In the final 
section of this paper we study how dissipation infiuences 
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spatial and temporal coherence in a finite size condensate 
and show how the linewidth of emission from polariton 
or atom condensate should be determined taking proper 
account of the spatial fluctuations. We further empha- 
sise the fundamental difference between emission from a 
polariton condensate or an atom laser and that from the 
photon laser. 

The paper is organised as follows: The model for the 
system, and for the reservoirs to which it is coupled is 
introduced in Sec. [Ill then in Sec. IIIII we show how to 
integrate out first the reservoirs, and then the fermionic 
fields to give an effective action in terms of the photon 
field. We then study this effective action in the saddle- 
point approximation in Sec. IIVI In Sec. |Vl by discussing 
fiuctuations about the saddle point we consider the sta- 
bility of the saddle-point solutions, and show how the 
instability of the normal state, and the photon distri- 
bution functions, compare to an equilibrium treatment. 
Having identified the stable and unstable saddle-point 
solutions, Sec. IVII then presents numerical results for the 
critical conditions at which steady-state, non-equilibrium 
condensation occurs. The effects of fluctuations on corre- 
lation functions in the condensed case are studied again 
in Sec. IVIIl where care is taken to correctly describe 
phase fluctuations in the broken symmetry system. Sec- 
tion |VlTT] then studies how flnite size modifles correlation 
functions, and the relation between the previous results 
and laser theory. Finally, section HXl summarises our re- 
sults. 



II. MODEL 

Our Hamiltonian is 

H = Hsys + -f^sys.bath + -f^bath, 

where. 



(1) 



a p 

+ ^EE(-9",pV'p&J.«a+H.c.) (2) 



describes two fermionic species ba and aa, interacting 
with bosonic modes ipp normalised in a 2D box of area 
L^, with L — *■ oo. Condensed solutions of Eq. ([2]) have 
been studied in the context of atomic Fermi gases^SiiSiii 
and microcavity polaritons^i^i^. In this work we fo- 
cus on microcavity polaritons, and so this model de- 
scribes the interaction between disorder-localised exci- 
tons which are dipole coupled to cavity photon modes 
■0p, with low p dispersion, tjp ~ + P^/2TOph, where 
Wph = {h/ c){2iT /w) is the photon mass in a 2D micro- 
cavity of width w. The disorder localised excitons are 
described here as in previous works^SiiSid^ by hard-core 
bosons; i.e. the Coulomb interaction between excitons is 
described by exclusion, preventing multiple occupation 



of a single disorder-localised state a. This hard core bo- 
son is represented by a two-level system, described here 
as two fermionic levels, hl^,a^. Thus, the combination 
fojjfl^ creates an exciton in the localised state with en- 
ergy Ca- This energy includes the Coulomb binding 
within an exciton state. In such a description, it is im- 
portant not to confuse the fcrmion states (representing a 
hard-core bound exciton) with the unde rlying conduction 
and valence band states (see e.g. Refs. I43ll45l for further 
discussion of this point). In order that these fermionic 
levels describe a two-level system, it is necessary that the 
constraint h\jD^ -\- a]^a^ = 1 is satsified; i.e. that exactly 
one of the two levels is occupied. In thermal equilibrium, 
this constraint can be exactly imposed by a shift of Mat- 
subara frequencies^; and in that case it can be easily 
seen that the difference between imposing the single oc- 
cupancy constraint exactly and imposing it on average 
leads only to a factor of 2 in the deflnition of tempera- 
ture. Out of thermal equilibrium, no simple shift to the 
Matsubara frequencies is possible, although an extension 
to the non-equilibrium case has been proposed^. For 
simplicity, in this work, we will impose the single occu- 
pancy constraint on average, as discussed below when 
introducing the occupation functions of the bath. 

Because of the imperfect reflectivity of the cavity mir- 
rors, photons escape, so the system must be pumped (ex- 
citons injected) to sustain a steady-state. As illustrated 
schematically in Fig. [1] the imperfect reflectivity of the 
mirrors is represented by coupling to the continuum of 
bulk photon modes. Incoherent fermionic pumping is de- 




FIG. 1: (Color online) Schematic diagram illustrating parts of 
the Fermi-Bose system, and its coupling to baths. The parts 
in the box labelled system are described by Eq. ((2)1 , while the 
effective couplings to the baths, described by Eq. ^ , lead to 
effective pump and decay rates 7 and k as discussed later. 



scribed by coupling to a pumping bath — which is repre- 
sented mathematically as two separate fermionic baths, 
coupled to the two fermionic modes. Thus, the coupling 
of the system to these pumping and decay baths is writ- 
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ten as: 

^sys,bath = E ra,fe(4 A- + H.C.) + r^,,(6ii3fc + H.C.) 

+ ECp.fc(^^*fc + H.c.), (3) 

p,k 

Here Af^^B^. are fermionic annihilation operators for the 
pump baths, while 5*^, are bosonic annihilation operators 
for photon modes outside the cavity. The Hamiltonian 
corresponding to the evolution of these baths is given by: 

iJbath = Y.^'i^ BiBf.+Y.^i^i'^f,. (4) 

The pumping bath, if thermalised at some finite non- 
zero temperature, acts both as a source of particles, and 
also tries to drive the polariton distribution function 
towards a thermal distribution in equilibrium with the 
bath. In some physical systems, one might also consider 
a bath which purely provides a thermalisation mecha- 
nism, such as phonons, which redistribute energy, but 
do not change particle number. We do not explicitly 
consider such a bath. However, in the example of micro- 
cavity polaritons, our model may still capture much of 
the important behaviour, for the following reason. One 
may consider the low energy polaritons as being pumped 
by a reservoir of higher energy excitons. These excitons 
are formed by the binding of the electrons and holes 
injected by the pump laser, and subsequent relaxation 
by phonon emission, and are thus partially thermalised. 
By regarding our pumping bath as describing a partially 
thermalised exciton reservoir, our model, being interact- 
ing, could thus describe the thermalisation of low energy 
polaritons pumped by such a reservoir. 

Although in the absence of other processes, the exci- 
tons would thermalise to the pumping bath, they are also 
strongly coupled to photons, which in turn couple to a 
second environment of the bulk photon modes outside 
the cavity. The strongly coupled exciton-photon system 
would be therefore influenced by two independent envi- 
ronments which are not in thermal or chemical equilib- 
rium with each other. Even in the steady-state, if the 
rates of dissipation to the environment are larger than 
the polariton-polariton interactions, the system would 
remain out of thermal equilibrium. In addition, even if 
the thermalisation via polariton-polariton interaction is 
fast, so the system distribution function would be close 
to thermal, particles are continuously added and removed 
from the system. We show that this particle "current" 
has dramatic consequences on the properties of such a 
condensate even if it remains close to equilibrium. 

We would like to stress that there are two distinct is- 
sues, both of which we intend to address. The first is 
that of non-equilibrium distribution functions, in systems 
where the internal thermalisation rate is slower than the 
pumping and decay rates — i.e. when the coupling to the 
external baths is strong, and the baths are not in equilib- 
rium with each other. The second issue is the presence of 



particle "current" in strongly dissipative systems — even 
if internal thermalisation rates are large, this current may 
be important if the pumping, decay and thermalisation 
rates are large compared to other energy scales. 

In the next section we will introduce the path integral 
formalism which will allow us to treat the nonequilibrium 
conditions. Our approach will then be to assume that the 
pumping and decay baths are much larger than the sys- 
tem, and so the populations in the baths are fixed. This 
will enable us to describe the properties of the system as 
influenced by its coupling to the baths. These influences 
modify both the system's spectrum and the population 
of this spectrum. We will look for steady states of the 
system in the presence of pumping and decay, and study 
the excitation spectra around these steady states. 



III. PATH INTEGRAL FORMULATION 

In order to study the system away from thermal equi- 
librium, we proceed using the path-integral formulation 
of non-equilibrium Keldysh field theory, as described in 
detail in Ref. |4^ Following the prescription there, we 
write the quantum partition function as a coherent state 
path integral over bosonic and fermionic fields defined on 
a closed-time-path contour, C. Arranging the fermionic 
fields into a Nambu vector = (6, a) and (j) — {h,a)^ , 
loosely referred to as "particle/hole" space, the partition 
function can be formally written as: 

p a 

k 

where M represents a constant of normalisation and the 
total action can be separated into constituent compo- 
nents S* = 5^ + 5.0 + S'batiK^ + Sha,t\i.iji- The part: 

S<P = J dt'^4'a [idt - daa^ ~ ga.p1pp(^- - .9a,pV'pO-+] 0a, 
C "'P 

describes the free exciton evolution together with the 
dipole interaction between excitons and photons. Due 
to the Nambu formalism, the term in brackets is a ma- 
trix, and has been decomposed in terms of the Pauli ma- 
trices (Ji operating in the particle-hole (6, a) space (with 
tTo = 1). The time derivative is taken along the Keldysh 
contour C. Similarly, 

^ dt '^{pp {idt - t^p) V'p 
c P 

describes the free photon dynamics. The excitonic envi- 
ronment and the interactions between excitons and their 
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environment is given by 



5, 



bath,0 



-Tl^j,{KBk + Bkb^) - K kiaaAk + Akac.)], 



while the photonic environment is given by 



5, 



bath,-?/; 



Cp,k (V^p^fe + ^'feV'p)]- 

As described in Ref. |48|, the standard procedure is to 
replace the fields on the closed-time-path contour by a 
doublet of fields tp = i/'b) on the forward and back- 
ward branches. This then leads to four Green's functions: 
forward iG<(t,t') = {iJf{t)i'l{t')), backward iG'>(t, i') = 

(V^6(*)V'/(i')), time-ordered iG'^{t,t') = {i, j{t)il}\{t')) , 

and anti-time-ordered iG'^{t,t') = I^i the 

homogeneous steady-state these arc functions of r — r' 
and t — t' alone and when transformed into p and lj 
space, in the case of photon fields, the functions iG"^ and 
iG^ give the luminescence and absorption spectra re- 
spectively. Again following Ref. lH, as these four Green's 
functions are not independent, one proceeds by making 
a rotation to classical 4>ci = {4>f + cj)b)/^/2 and quantum 
<j)q = {(j)f — (j)},)/\/2 components. All fields are from now 
vectors in Keldysh space, i.e if) = {tpd, ''Pq)j ^.nd we define 
an additional matrix in Keldysh (d, q) space: 



M 



1 f ^cl -09 

V2 V V'c/ 



(where are Pauli matrices in Keldysh space). One 
may then write the action as: 



dt 



<Pa [idt 



a,p 



9a,p I M 1 K i 



/oo 
di X! ^ ^1 ^P' 

"°° p 

/oo ^ 

-°° a,k 
- ^a,k{daCf\Ak + AkO^aa) 

+ Bk {idt - 4") <ylBk + Ak (idt - 4°) <^IM, 

/OO 
dt^[-Cp,k (V^pcr^^fc + ^kcrilfjp) 
-OO __ 



p,A: 



A. Treatment of environment 

As we are interested in the properties of the system, 
rather than the properties of the baths, we next inte- 
grate over the bath fields, to leave an effective action 
expressed only in terms of the fields describing the sys- 
tem. If the baths are much larger than the system, then 
their behaviour is not affected by the interaction with the 
system. One may then evaluate correlation functions of 
bath operators as for free bosons and free fermions; these 
correlators in turn depend on the distribution function of 
the baths, i.e. the population of the bath modes. The 
effects of the environment then enter as self energies for 
the system fields, which modify both the spectrum and 
its occupation. This procedure is described in Ref.lH; we 
summarise the results here both to show how it applies 
to our system, and also as our notation differs slightly 
from Ref. IH For the decay (photon) bath one has: 



Shathip — / dtdt' ^ 

■J-^^ p.p' 
MthiT.^p^kCp'M [iidt-4y( 



'<jf^p'it'), 



In Keldysh space the Green's function for a free boson 
has the following form 



[idt 



(t-f) 







where (after the Fourier transform with respect tot — t') 
the retarded, advanced and Keldysh Green's functions 
are respectively 



D 



RjA 



1 



If the bath distributions are thermal, then tib would be 
the Bose occupation functions, however one can also con- 
sider arbitrary function for ns- 

Let us now make a number of restrictions on the pho- 
ton bath, to simplify the analysis. Firstly, we will as- 
sume that 5'bathV' does not contain terms off-diagonal in 
pp'. This means that each confined photon mode p cou- 
ples to a separate set of bulk photon modes, i.e. that 
Cp,fcCp',fc = unless p = p'. Physically, this can be in- 
terpreted as conservation of in-plane momentum in the 
coupling of two-dimensional microcavity photon modes 
to bulk modes. Next, we restrict to the case that all 
p photonic modes couple to their environments with the 
same strength i.e Cp^^ = C^. Then, if the bath frequencies 
cij^ form a dense spectrum, and the coupling constants Cfc 
are smooth functions of the frequencies, we may replace 
the sum over bath modes by an integral, 



fc 



duj'^Ci^'^fN'^i'^^), 
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where we have introduced iV^(a;'') as the bath's density 
of states. After integrating over uj'' we obtain^ 



'bath'i/' — 



ipp{-uj). 



i-co) 



By writing d-^'"^{uj) = R{Lu)^iK{ijj) we may spHt the bath 
self energy into an imaginary part, describing broadening 

and a real energy shift, 



Riuj) = J d^''- 



In terms of these, the Keldysh component becomes: 
d^{u}) = -i2K{ui){2nB{uj) + 1). 

Although the formalism allows one to consider any 
density of states, and coupling strength as a function of 
frequency, one possible choice is a Markovian (or Ohmic) 
bath — i.e. a white noise environment — where the den- 
sity of states for the bath and the coupling constant of 
the system to the bath are frequency independent, and 
so C^(a;^)-/V^(a;'') = ('^N'^. For this case the real en- 
ergy shift R{io) is zero while k{u!) = k. In this work, 
we will consider this Markovian limit, but due to the 
bath's occupation function, the Keldysh component will 
remain frequency dependent. Combining the free pho- 
ton action with the effective action for the photon decay, 
using F^(aj) = 2nB{ijj) + 1, one has: 



/oo ^ 
dui 
-oo „ 







— iO-n 



■0p(-Cj). 



One can follow a similar procedure for the baths con- 
nected with the pumping process. 



5, 



bath0 



f f dtdt'Y^ 



k 
k 

The Green's function for a free fermion is 



P^{t~t') 



where in frequency space 



1 



In the same way as above, np would be the Fermi occu- 
pation function for a thermal distribution. 

For compact notation, we will define additional matri- 

... /^io\, /OO 

CCS m (b,a) space as CT| = I q g I and (T| = 



and so: 



1 



E 







where n^/^/^ - V r^'^^r''/" with the 

fermionic propagators p^^^/^ as defined earlier. Now 
we make similar restrictions as for the photonic environ- 
ment: we consider all excitons coupled equally strongly 
to the environment (the coupling constants of the system 
to the bath is a independent) and take the Markovian 
limit. Without much loss of generality we can further 
assume that the coupling strength of the two fermionic 
species to their pumping baths are the same, after all of 



which r 



b/a 



r. As in the bosonic case, in the Markovian 



limit, the real self energy shift vanishes, and imaginary 
part takes the form: 
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ttF^TV^ 



with Np being the bath's density of states. Of course, 
due to the distribution function of the bath, despite the 
Markovian limit, the effective pumping rate of a given 
exciton state will depend on its energy. The final form is 
then: 



-i7cro 
^70-0 2i^{Fi,{~v)a^ + Fa{-v)(y{) 



where Fi,{v) = 1 - 2-n}'p{v) and Fa{v) = 1 - 2n'],{v) are 
the fermion distribution functions. 

Any functions for the bath distribution and can 
be considered within this formalism. One physical choice, 
as illustrated in Fig. [21 can be pumping of quantum-well 
excitons by contact with some thermal reservoir with a 
chemical potential /i^, i.e 

Fi,{v) = taiih^{v - ^ib) Fa{v) ^ t&n}i^{v + ^ib), (5) 

where (3 = 1/kT. Note that, as discussed earlier, these 
bath distributions have been chosen so that on average, 
{U^h + a) = 1; i.e. the single occupancy constraint for 
these fermionic states to represent two-level systems is 
obeyed on average. In the absence of any other processes, 
contact between the excitons and the pumping reservoir 
would control the population of excitons, and so 



(fetfo-a^a) = n^(e) - n^(-e) 



-tanh 



2 



Mb)- 



Thus, by pumping with a thermalised source of electrons, 
one will find a thermalised distribution of excitons. 
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FIG. 2: (Color online) Schematic diagram illustrating occupa- 
tion of A and B baths, and their coupling to the two fermionic 
levels. In the case shown, the bath chemical potential is below 
the energy level e, so the pumping cannot lead to inversion. 
The flow of energy is from the pumping bath, through the 
fermionic levels of an exciton, to the photons, and energy is 
then lost into the photon bath. 



ciple one could consider a non-Markovian environment 
where for some range of frequencies one has 7(0;) = 0, 
k(cj) = 0. Then, for that range, there would be no 
damping, but also (lu) = 0, {uj) = 0, so the sys- 
tem distribution in such a range would not be influenced 
by the bath — i.e. no thcrmalisation. Thus for a full 
thcrmalisation of all relevant modes of the system, one 
needs a non-zero coupling to the low frequency modes of 
the environment, which will at the same time introduce 
dephasing. 



B. Integration over fermionic fields 

After eliminating the bath's degrees of freedom the full 
action S becomes 



S 



dtdt' 



Before proceeding further, let us examine what the 
form of the self energy due to the bath tells us about the 
relation between thcrmalisation and dephasing. In prin- 



where G is the exciton Green's function. Introducing the 
abbreviations A^; = ^'^p.d ^nd A, J2p ^'^p.i 
we may write G as: 



{Xq{t)a^ + Xq{t)a_)Sa,a' 



- Ac;(i)CT_|_ - Xcl{t)a_)Sa,a' - ijCTa 



Acz(t)cr+ - Xcl{t)cr^)Sa.a' + ijCT^ -(Aq(t)cr_|_ + \q{t)(J_)5 a.a' + 2i7(FbCr^ -j- FaCJ ^) 



(6) 



Note that Fbit — t'), Fa{t — t') as well as idt are time 
non-local. It is now explicit how the competition be- 
tween the two environments works. The photon environ- 
ment, with the distribution function F^ of modes outside 
of the cavity, affects the free photon evolution. Simi- 
larly, the fermionic environment, with the bath distri- 
butions Ff),Fa, enters the exciton Green's function, now 
also modified by the presence of cavity photons. The 
spectrum of this coupled system will combine both the 
strong coupling between excitons and photons as well 
as the dissipation to the environment. The occupation 
of these modes will be a non-trivial combination of the 
distributions of the baths as well as the exciton-photon 
interaction. 

The action S is quadratic in fermionic fields and there- 
fore it is possible to integrate over the fermionic degrees 
of freedom (jj and, in the same spirit as in previous stud- 
ies of the equilibrium properties of this model^^i^i^i^. 



obtain the total effective action for the photon field alone 



^ = -i^TrlnG^i„+ / dtdt'Y^ 



Other than those approximations explicitly discussed in 
the text, this expression is exact; i.e. it makes no assump- 
tion about what form '0p(i) takes. Note, however, that 
due to the non- linear term TrlnG|J^^ the action is highly 
complex and contains all powers of ^ci and ipq ■ Therefore 
some expansion scheme needs to be performed. 
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IV. SADDLE-POINT (MEAN-FIELD) ANALYSIS 

In order to determine the state of the pumped, decay- 
ing, strongly coupled system, we will follow a standard 
method for path integrals, and first find the saddle point 
solution. The saddle-point equations for the action ([7]) 
have the following form 



6S 



= dt'\ [{idt - u;p)S{t - t') - cf {t - t')] Vv,c/(i') 



-d^(t-OV'p,,(i')}-E ^(-*)Tr(G<,,„a_ao^) = 0, 



5S 

Si>p,ci 



dt' [{idt - UJp)d{t - 1') - d^{t - 1')] ^p,,(t') 



-E7fH)Tr(G„,„a_af) = 0. 

It can be seen that the second equation is satisfied by 
V'p.ij = (classical saddle-point). Putting i^p^q = into 
([S]) gives the usual structure for the mean-field exciton 
Green's function, which ensures causality^: 



and so G 



G-1- 



G^ G« 







G^]-i 

K 



, where G^ = -G"[G-^]'^G' 



G^ 

It is now clear why the Keldysh rotation discussed ear- 
lier, i.e. working in terms of (d, q) components rather 
than {f,b), is more convenient. By reducing the number 
of dependent functions, both the Green's function and 
inverse Green's function contain a zero block, and so be- 
come easier to invert. With this structure Tr(G(7_(To) 
Tr(G^cr_) and Tr(Ga_a5^) = Tr((G^ + G'^)a_) = 
since G^{t,t) + G'^it^t) = 0. Thus we are left with only 
the first of the saddle point equations, which now be- 
comes: 



dt' [{idt - ujp)5t-t- - d'^it - t')] ijpAt') = 

5:^H)Tr(GiVj. 

Since we consider an infinite homogeneous system (no 
trap) we expect a uniform saddle point. We therefore 
consider the solutions to be of the form ijjp = ip6{p). It is 
difficult to invert [G~^]^'^ matrix for an arbitrary time 
dependence of the ip fields. We are, however, interested 
in the non-equilibrium steady-state so we take the only 
time dependence of the photon field ip to be oscillation 
at a single frequency. Therefore we propose the following 
ansatz 



Substituting this ansatz in the action of Eq. ^ will 
lead to explicit time dependence within the exciton in- 
verse Green's function. This time dependence can be re- 
moved straightforwardly by implementing an appropriate 
gauge transformation, described by the following matrix 
in particle-hole space: 



U = 











The trace is invariant under unitary transformations, 
TrlnG^i = TrlnUG'^W , so the efl?ects of such time- 
dependence appear in only two places: Firstly in the 
time derivative terms, which lead to the energy shifts 
LOp ~* ujp = ujp - ^J,s and ^ = e« - ^J■s/'2, 
and secondly in a gauge transformation of the bath 

functions rf«M/^(u;) ^ rf«M/^(u; + fj.s), V^^^^^iy) ^ 
pfM/^(,+ i^) and pfM/^(^)^pfM/^(.-i^).In 
practice, the latter substitutions mean replacing I-Ib ^ 
fj-B ^ fJ-B - l^s/^ in Fa, Fb. 

With tpq = and the above time-dependence described 
by the gauge transformation, the matrix G~^ can now 
be easily inverted and the final form for the mean-field 
exciton Green's functions are 

i^^ — E'^ ± 2i"fi> — 7^ ' 
Fb/a{i^)[{,.±~ea?W] + Fa/bii^)g'm' 



and 



[(z.-i?„)2+72][(^ + iJ^)2+^2 



GtH = -(Gffc(i^))* = -2j7#/x 
where a, b defines the particle-hole space as follows 

Gbb Gba 



(10) 



(11) 



G = 



Gab Ga 



and Ea ~ \/^a + S^IV'/P- Note that only the site-index 
diagonal, i.e. (a, a), component appears in the gap equa- 
tion and so we have omitted the site index in G for 
brevity. Also since at the saddle point t/jq ~ then 

V2 



In this work we consider gp^a = 9- 



The influence of the distribution of the oscillator strength 
has been addressed in Ref. |4^. The mean-field exci- 
ton Green's functions physically correspond to excitons 
strongly renormalised by the presence of the mean-field 
photon field, and damped by the coupling to the envi- 
ronment. The Keldysh Green's function, which contains 
the distribution of excitons, depends on the distributions 
of the pumping bath. In general 



^(t) = V^e"*^^*. 



(8) 



G' 



G"F - FG^, 
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where F has a meaning of the quasi-particle distribution 
function. We can determine the mean-field distribution 
function for excitons in a self-consistent photon field from 
Eqs. ®-(in]): 



bb/aa 



(^) 



F,{v) ^ {F,{u) 



Fail^Wl+j') 



FM = {FabTil^) = 



g^f{Fb{v)~Fa{y)){i^+i-f) 



where Fa and F^ are the bath's distributions given by 
Eq. (O, with ^ A'-B- Note that, since this is a mean- 
field approximation, only coherent photons enter in this 
distribution. Thus, in the uncondensed case where V' = 
the cxciton distributions reduce to -Fh(,/aa = ^fc/a and 
Fab = 0: i.e. in the absence of coherent photons, the 
mean field approximation neglects the effect of photons 
on the exciton distribution. The distribution of the pho- 
tonic environment will however enter the distribution of 
fluctuations about the mean-field, as will be discussed in 
SecES 

With the ansatz ([5]) the saddle-point equation becomes 



(wo - MS - iK)i'f 



y-' 

^ 2 



(-z)Tr«J. (12) 



As in equilibrium this is a self-consistent equation for 
the order parameter (condensate). With changing den- 
sity the type of transition moves from interaction domi- 
nated BCS-like mean-field regime to a fluctuation dom- 
inated BEC hmit^ (strictly speaking BKT in 2D) . So 
the above equation is analogous to the gap equation in 
the theory of BCS-BEC crossover. Here it relates the 
coherent photon field with the exciton Green's function 
strongly modified by the presence of such a coherent field. 
Physically it means that the coherent field is generated 
by a coherent polarisation in the exciton system which 
in turn is generated by the presence of the coherent field. 
Thus equation \Y1\ can be viewed as a non-equilibrium 
generalisation of the gap equation. One difference with 
respect to equilibrium is that the distribution function 
contained in now may not be thermal. However, the 
more important difference is that the gap equation (|12p is 
now complex and gives two equations for two unknowns: 
the order parameter '\\) and the frequency \is- 

The common oscillation frequency \xs would in ther- 
mal equilibrium be the system's chemical potential, con- 
sidered as a control parameter, adjusted to match the 
required density, and the (real) gap equation determines 
only ?/'. Here, because different baths have different 
chemical potentials, the system is not in chemical equilib- 
rium with either bath, so both ^5 and i/" must be found 
from the gap equation. The density, which can be found 
given ?/) and /is, is set by the relative strength of the 
pump and decay. 

Thus the real part of the gap equation is analogous 
to the gap equation for closed equilibrium system, where 



the right hand side describes polarisation due to nonlin- 
ear susceptibility. By considering the existence of pump- 
ing and decay, one also introduces the imaginary part, 
which describes how the gain balances the decay (as in 
lasers) but now in the strongly coupled exciton-photon 
system. If one were to instead consider the equilibrium 
theory, and merely add decay rates, one could not a priori 
guarantee that the fluctuation spectrum would be gap- 
less, as should arise from spontaneous symmetry break- 
ing. By ensuring that gain and decay balance, the fluc- 
tuation spectrum above the ground state which satisfies 
both the real and imaginary parts of the gap equation, 
will indeed be gapless. By connecting the equilibrium 
self consistency condition (gap equation. Gross Pitacvskii 
equation), and the laser rate equation, Eq. (fT^ puts the 
condensate and the laser— in the same framework and 
so allows study of the crossover and the relation between 
the two. 

Using pT|) the mean-field equation becomes 

(wo - A*s - jK)'0/ = 



2 (Fa+fb)^+(Fb-fa)(ga+Z7) 

27r ^ ^ [(:/ - i?„)2 + 72] {(y + E^f+ 72] • ^ 



Note that 7 appears both in the denominator (as it gives 
rise to the dephasing) and in the numerator (it gives rise 
to pumping). 

As in thermal equilibrium, the normal state V'/ = 
is always a solution of Eq. (|13p . but for some range of 
parameters there is also a condensed ?/;/ 7^ solution. 
For ipf ^ the final form of the gap equation is 



UJQ — IK 



9^1 



E 



du {Fa + Fb)iy + {Ft - Fa){i^ + Z7) 
27r [(:/-i?„)2+72][(,. + £;„)2+72] 



(14) 



Now for a given set of parameters we can solve the real 
and imaginary parts of this equation to determine the 
coherent photon field ipf and its oscillation frequency /is- 

We can reduce the number of parameters in our the- 
ory by measuring energies in units of the exciton-photon 
coupling g and, noting that our equations have made 
no assumption about the origin of energies, taking eg 
as some reference energy, such as the bottom of the ex- 
citon band. The independent parameters in our the- 
ory arc then the distribution of exciton energies [i.e. 

—^ J dei>si^)]j the detuning of the photon from the 
reference point A = ujQ — 2eo, the pumping bath chemical 
potential /is — eo, the pumping (decoherence) strength 7, 
and the coupling to the decay bath n. 

Having found the self consistent oscillation frequency 
fis and coherent field, one can then calculate the ex- 
citonic density and polarisation. The polarisation, i.e. 
(a'b) (where |(a''^6)|2 also gives the number of condensed 
fermion pairs - condensed excitons) follows directly from 
the gap equation, so the magnitude of polarisation is 
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given by \JilP- ^ n'^ipf. The excitonic density is given by: 

a a 

(Fb ~ Fa) [g^lV^p - - il - 7^] - (ffc + Fa)2i,i^ 

(15) 

Since our choice of bath populations in Eq. ([5]) imphes 
that the empty state corresponds to {b^b) = 0, {a^a) = 1, 
it will be convenient to shift the exciton density so that 
the empty state corresponds to zero density, thus: 

Poxciton '^'Y^{^+ b^ba - a^aa) ■ (16) 
a 

In the limit that the temperature of the pumping bath 
goes to zero, one can perform the various integrals in 
Eq. (|14|1 and Eq. (|15p in terms of elementary functions. 
These forms are presented in Appendix [Al 



Let us consider the limit where k/^ —f 0, i.e. coupling to 
the photon bath vanishes faster, and so the distribution 
will be set by the pumping bath. Then the left hand side 
of Eq. pB]) is zero and so one requires Fa{E) = Fb{E). 
Using the gauge transformed versions of the thermal dis- 
tribution functions in Eq. ([S]), this condition becomes 
jls ~ 0, i.e. that fj,s = 2fj,B- Putting this solution into 
P7p we recover the equilibrium gap equation at a tem- 
perature T set by the pumping bath: 

.9' P 
^0 = — tanh-E. 

This limit provides a reassuring test of the formalism, 
and also supports the interpretation that the real part of 
the gap equation connects the order parameter with non- 
linear susceptibility, while the imaginary part describes 
the balance of gain and decay, and so controls /ig and 
the particle density in the system. 



V. SECOND ORDER FLUCTUATIONS AND 
STABILITY OF SOLUTIONS 



A. 7 = 0, K = limit 

In order to understand the meaning of the gap equa- 
tion, and the connection to condensation in a closed 
equilibrium system, it is instructive to take the limit 
7 ^ 0, K — > in Eq. p^ . This will also provide a con- 
sistency check of the non-equilibrium theory as it should 
recover the equilibrium limit as the coupling to the envi- 
ronment approaches zero. The real part of Eq. can 
be rewritten as: 



Wo - Ms 



9 

AE 



dv 
2^ 



7 



7 



{v-EY+^'^ {v + Ef+^'^ 



{Fa{y) + Fb{v)) + (Fbiiy) - Fa{v))^ 



From the definition of the 5 function we have 
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lim — .„ „ 

7-0 {v - EY +7^ 



tt6{v-E) 



and so, using Fa{—v) = —Fb{i'), the real part of the gap 
equation reduces to 



UO - M5 



AE 



[FbiE) + Fa{E)] 



^[Fb{E)~Fa{E)]. (17) 



Similarly, the imaginary part of the gap equation can be 
rearranged as: 



K 

7 



4E^ 



[Fa{E) ~ FbiE)] 



(18) 



Having found the self-consistency condition, consider- 
ing the possibility of uniform condensed solutions, we 
next consider the stability of such solutions. The con- 
sideration of stability is important firstly since, as dis- 
cussed above, ipf = is always a solution of Eq. ([T^ . so 
one must determine which of the normal and condensed 
solutions is stable, and secondly because we considered 
only spatially homogeneous fields, with a single oscilla- 
tion frequency, so one may find that neither ipf = nor 
our ansatz of Eq. ^ is stable, suggesting more inter- 
esting behaviour. There is an important difference in 
interpretation of the saddle point equation between the 
closed-time-path path-integral formalism used here, and 
the imaginary-time path-integral in thermal equilibrium. 
In the imaginary time formalism, extremising the action 
corresponds to finding configurations which extremise the 
free energy; thus, stable solutions correspond to a mini- 
mum of free energy, and unstable to local maxima. Here 
in contrast, for a classical saddle point (i.e. ipq = 0), 
the action is always S = 0, and the saddle point condi- 
tion corresponds to configurations for which nearby paths 
add in phase. Thus, in order to study stability one must 
directly investigate fluctuations about our ansatz, and 
determine whether such fluctuations grow or decay. 

In considering the question of stability, we will first dis- 
cuss stability of the normal state, which is instructive as 
it shows how the question of whether fluctuations about 
the non-equilibrium steady-state grow or decay is directly 
related to the instability expected in thermal equilibrium 
systems when the chemical potential goes above a bosonic 
mode. We will then turn to the spectrum of fluctua- 
tions about our condensed ansatz. While we will discuss 
here whether such fluctuations are stable or unstable, 
we will defer until Sec. IVIII the evaluation of correlation 
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functions associated with these fluctuations. This is be- 
cause, as discussed there, these fluctuations include phase 
modes, and phase fluctuations may become large. It is 
therefore insufRcient to only expand to second order in 
fluctuation fields, but one must instead reparameterise 
ip = y/pQ + 7re"^, and then describe the correlation func- 
tions of tjj in terms of those of phase and amplitude tt, 
including the effects of 4> to all orders. Such a compli- 
cation is however not needed in order to study whether 
fluctuations are stable or not, and so it is reasonable to 
postpone such a treatment, and consider an expansion in 
terms oi ip = ipo + 6ip to second order in Sip instead. 

Thus, to find the spectrum of fiuctuations, we consider 
the effective action governing fluctuations about either 
-0 = V'o or about ■0 = 0. Considering the effective ac- 
tion in Eq. ([7]), and expanding to second order in Sip, 
one finds a contribution from the effective photon action, 
and a contribution from expanding the trace over exci- 
tons. This latter contribution can be found by writing 
Ga!a = (G^'^a)"^ + SGa \ where G"P is the saddle point 
fermionic Green's function, which depends on the value 
of the saddle point field ipf, as given in Eq. (fTU]) 
and mi), and the contribution of fluctuations 6G~^ is 
given by: 



-^{gSipcia- 



Thus, one can expand the action as: 



- ^ Trln [G-^ = (-^) ^ Trln [(G^p J-^] 

a. 

+ (-*)(-^)ETr [GZJG-^'GZJG-'] . 

a 

In this expansion, we have retained only the terms diag- 
onal in site index; i.e. neglected any bath induced inter- 
action between different exciton sites. Such bath induced 
interactions should be small for small 7, and their inclu- 
sion would considerably complicate the formalism. Such 



an approach is also equivalent to considering a separate 
set of baths for each disorder locahsed state a. 

Because, in the presence of a coherent field, the ef- 
fective action can contain terms like SipS'tp and SipSip, it 
is convenient to introduce a Nambu structure of photon 
fields. Thus, the photon fiuctuations are described by a 
2x2 = 4 component vector, with one factor of 2 from the 
Kcldysh structure, and one from the Nambu structure, 
hence: 



SA = 



( 6ipci{uj) 

\ 5%l)q\~Uj) 

in terms of which the action for fiuctuations 5Sf is: 



(19) 



5Sf 



2^ 



V' 



[V- 



1 K 



5K(uj). 



For convenience later, we shall introduce the notation: 



_11 R/A/K 



R/A/K j.B,/A/K 



K. 



R/A/K j^R./A/K 



(20) 



By definition we have that: \D-^\^ = ([l?"^]'") , and 

in addition the Nambu structure implies certain symme- 
tries between the elements of \D~^\^^^I^ , which together 
can be written as: 



i^f(c.)=i^f(-c.) 



Introducing the compact notation: 
dv 



(21) 



(22) 



[1*9 
we may thus write: 



2tt 



f{v)g{v~uj) 



' —oj — w, 



■ 9 



'-'bb * ^aa + ^bb * ^aa '-'ba * ^ba + ^ba * ^ba 



ab * ^ab + ^ab * <-^a6 



'^aa * '^bb 



^aa * '^bb 



(23) 
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and 



2iKF,f,{—uj + fis) 




ab 



ab 



ah 



(24) 



r<A ^ c'R 

'^aa * '^bb 



Normal state excitation spectra and 
distributions 



Thus, the normal state luminescence, absorption, and 
distribution functions are given by: 



The excitation spectrum can be found from the poles 
of the fluctuation Green's function, i.e. from the zeros of 
det To extract the occupation of the spectrum, 

one can extract the boson distribution function via 



^ _^R 



where simply T>^l^ = [l? ^] . Whilst in general 

these are 2x2 matrices in Nambu space, in the normal 
state this structure is redundant, and so the distribution 
fimction is the diagonal constant matrix Fg = 271^ + 1, 
where ns describes the occupation of the modes. Alter- 
natively, one can invert the Keldysh rotation in order to 
find the physical Green's functions. 



F.V 



(25) 



which as discussed in Sec. IIIII relate directly to the lu- 
minescence, £(tj,p) = il?<(a;, p)/27r, and absorption 
^(WjP) = i2?>(a;, p)/27r. Still, assuming the normal 
state, so that the Nambu structure is redundant, these 
become: 

I?«(c.,p) 



£(a;,p) = ns(a;)Im - 



A[uJ,^) = (ns(a>) + l)Im 



P«(co,p) 



While this form illustrates how the spectral weight and 
occupation can be separately extracted from the lumines- 
cence and absorption, in order to study these quantities 
it is more helpful to write them in terms of the compo- 

R FC 

nents, A'j^ ' of the inverse Green's function. In the nor- 
mal state, there are no anomalous (off diagonal in Nambu 



space) contributions, and so /Vj 



R/A/K 



K. 



R/A/K 



0. 



{C,A) (c^,p) = 



Fs{uj) 



-iK^{uj) =F 2Im [K?iuj) 
47r|iff(^,p)|2 

-iK^{u) 



(26) 



(27) 



2Im [Kf (w)] ' 

Let us now discuss what can be understood in general 
from the form of these equations, and then illustrate this 
discussion with the simple case 7 ^ T. From the differ- 
ence of luminescence and absorption in Eq. (|26p. one can 
identify a spectral weight: 



27rS'(tJ,p) = 



Im [A^f (w)] 



Rc[/vf(w,p)]'+Im [K^{uj)Y 



(28) 



thus, if the imaginary part of K^'{lo) is a smooth function 
of omega, then one will have almost Lorentzian peaks of 
the spectral weight at values uj* where Re [isrf-((jj*)] = 0. 
The width of these peaks, i.e. the lincwidth, is then 
given by Im [ii'f^(w*)]. Thus, the imaginary part plays 
one role as determining the linewidth. It also plays a 
second role, since from Eq. (|27p . a zero of the imaginary 
part causes the distribution to diverge; however, at these 
same points Eq. ([2H1) implies the spectral weight vanishes, 
so the number of photons does not diverge. Since a Bose 
distribution would diverge at the chemical potential, we 
can use this as a definition of an effective boson chemical 
potential, so Im [A'f^(/i'^^)] = 0. These results are illus- 
trated in Fig. [21 which show the luminescence, absorp- 
tion, spectral weight, and distribution function against 
the real and imaginary parts of and K'^ . 

From the above, it is clear that the form of Im [/vf^(ti;)] 
as well as K{^{uj) conspire to set the effective photon 
distribution. Using the expressions in Eq. (|10p . 
and (UJ), and for the moment restricting to the case 
= e we may write: 
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iK^{uj) = K,F^,{u)) + 



2Im [iff (w) 



I 2 2 



4 



2Rc 



dv 



2t: {v — e + — lu + e — 27) 



AY 



2-K [{v - e)2 + 72] - + g)2 + -y2] 

I — 



2^ [(i. - e~)2 + Y] + e)2 + 72 

(29) 



(30) 



For the case of pumping baths being individually in ther- 
mal equilibrium, one may get some insight into how the 
pump and decay baths compete to set the systems distri- 
bution. In the limit 7 ^ T, where T is the temperature of 
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FIG. 3: (Color online) Relation between real and imagi- 
nary parts of and and the luminescence, absorp- 
tion, spectral weight and distribution function. Plotted for 
7 = 0.2j(,K = 0.02g,T = 0.1ff,Ms = -O.Sj;. (cf Fig. H]) The 
upper panel shows the real and imaginary parts of J^f , and 
their zeros marked by arrows. The lower panel shows how 
these lead to peaks in the luminescence spectrum, and how 
the zero of Im[Jv^] defines divergence of the distribution func- 
tion. 

the pumping bath, the distribution functions Ff,/a{i^) arc 
smooth, while the denominators lead to sharp peaks, of 
width 7. One can then approximate the integrals by as- 
suming that over each Lorentzian peak, the distribution 
function takes its value at the maximum of that peak, 
and so: 



kFip (tj) 



Fsito) 



K + 



g^-f{l-Ftii)F4i-io)) 
juj- 2g)2 + 472 
g^^{F,{i)^Fa{i-uj)) 
{ui - 2?)2 + 472 



(31) 



From this one can see immediately two trivial limits. 
If 7 = or if g = 0, then there is no influence of the 
pumping bath and so Fs{uj) = Fq,(Lu), i.c the photon 
distribution in the system is the same as the distribution 
of bulk modes outside the cavity. Similarly, if k = 0, the 



photon bath has no effect, and 

1 - Fb{e)Fa{e - uj) 



Fs{uj) 



Fb{i) - Fail - LO) 



Thus, as one might expect, if the fcrmions arc in ther- 
mal equilibrium with Fb.a{i^) = F{v^ /is) where F{v) = 
tanh(/3i^/2), then by using a standard hyperbolic trigono- 
metric identity, this gives a thermal Bose distribution for 
the photons, with the same temperature, but twice the 
chemical potential, as expected since one boson corre- 
sponds to two fcrmions: 



Fs[ijj) = coth 
= coth 



- fJ-B) - ^ - UJ + ^IB) 



{u - 2//_b) 



The above expressions have been written after the gauge 
transformation described following Eq. Q. Of course, in 
the normal state, such a gauge transform has no effect, 
since it just corresponds to an arbitrary shift of the ori- 
gin for measuring energies, but we use the transformed 
notation for consistency with the condensed case. 

More generally, the two distributions compete to con- 
trol the photon distribution, which in general will not 
be thermal even if the baths are individually thermal, 
because they have different chemical potentials and tem- 
peratures. It is clear from Eq. ([?!]) that the effect of 
the pumping bath is largest near uj = 2e, and far from 
this value, both numerator and denominator are instead 
dominated by the photon bath. Physically, this means 
that the effect of the pumping bath is only important at 
energies where the photons are nearly resonant with, and 
so couple strongly to, the excitons. 



B. Instability of the normal state above the 
transition 

The discussion in the previous section, which defined 
H'^^ by zeros of the imaginary part of iff, and tj* by 
zeros of the real part allows one to understand the insta- 
bility of the normal state. It can be seen that the gap 
equation, ([T^ . if evaluated at "0/ = 0, is equivalent to 
the condition K^iuj = 0,p = 0) = 0, (measuring w rela- 
tive to fjLs)- This can be understood physically by seeing 
that the vanishing of A'f (w = 0, p = 0) implies there is 
a zero mode, corresponding to global phase rotations, as 
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one expects in a broken symmetry system. Thus, this 
condition imphes that there is a frequency at which both 
real and imaginary parts simultaneously vanish; i.e. the 



gap equation is the condition that ^' 



off 



the effec- 



tive chemical potential reaches the bottom of the normal 
mode spectrum. One can say "bottom of the spectrum" 
since the p dependence only enters the real part of K^, 
and Lo* will increase as p increases, thus if Wp^g < M''*^: 
then there will be a non-zero p for which cj* = ji"^ . Thus, 
the existence of a non-trivial solution to the gap equation 
can still be understood as a "chemical potential" reaching 
the bottom of the band^, even in this non-equilibrium 
context, as is illustrated in Fig. [H 
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FIG. 4: (Color online) Energy of zeros of the real and imagi- 
nary part of as a function of /is, demonstrating how 
the zero of Im(7('f) approaches oj*, a zero of Re [iff ] at the 
transition. Results are plotted for hi — 0.02g,TB = O.lg, and 
three values of 7 as indicated. The dotted vertical line marks 
the locations of the trace plotted in Fig. |3] 

It is also possible to connect the effective chemical po- 
tential reaching the bottom of the band to instability of 
the normal state, i.e. fluctuations growing in time. Let 
us consider poles, of the retarded Green's function, i.e. 
zeros of Ki {£^p,p). If these poles have negative imagi- 
nary parts they correspond to fluctuations that decay in 
time, and if positive, to growing fluctuations; thus sta- 
bility requires the imaginary part to be always negative. 
It is clear that at large enough momenta, the Green's 
function is that of bare photons, and is stable. Thus, 
if there are to be unstable modes, then there must be 
some p value at which the imaginary part of the poles 
goes from negative to positive. For reasonable systems, 
where the linewidth is a smooth function of momentum, 
this means the imaginary part must go through zero. A 
zero of Im[^p] means there is a real frequency which sat- 
isfies K^{£^p,p) = 0. However, the existence of a real 
frequency satisfying this condition was, as discussed pre- 
viously, exactly the gap equation at '0 = 0. Thus, if 
— oj^ — for some |p| = Pc, then for |p| < pc 
one will find positive imaginary parts. To illustrate this. 



consider a linear expansion in cj, so that: 

Iff (c., p) ~ (c^ - cj;) + ia{uj - M^^) ^ C{uj - Cp), 

then one finds, Im^p oc [fjp^ — cj*). 

Two more important connections can be drawn from 
the relation between poles of the retarded Green's func- 
tion, the distribution, and the gap equation. The first is 
that, as for any second-order phase transition, approach- 
ing the phase transition from the normal side, the fiuc- 
tuation Green's function describes a susceptibility which 
diverges at the transition. The second relates to the dual 
role that Im [A'f (w*)] played as the linewidth. As one 
approaches the phase boundary, at which real and imag- 
inary parts both have zeros, one must have that the ef- 
fective linewidth vanishes. These points are illustrated 
in Fig.O Note however that Im [i^/^(<x')] is of course not 
a constant, and so there will be some non-trivial line- 
shape, but a linewidth defined by full width half maxi- 
mum will vanish on approaching the condensed state, as 
a peak develops at = 0. The vanishing of homogeneous 
linewidth at the transition is a manifestation of diverg- 
ing susceptibility in an infinite system. Finite system 
size is expected to smear out this divergence and result 
in the homogeneous linewidth remaining non-zero, but 
still having a minimum near the transition. Additionally 
inhomogeneous broadening of exciton energies will add 
to the linewidth measured in experiments. 
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FIG. 



5: (Color online) Linewidth of the lower mode (solid, left 
axis), and energy integrated luminescence py = (dashed, 
right axis) in the normal state as a function of pumping bath 
chemical potential, as the phase boundary is approached. Re- 
sults are shown for two different dephasing parameters 7 with 
K = 0.02(7 for all three. 



C. Fluctuations in condensed state - stability and 
collective modes 

From the previous section we conclude that when there 
is a non-trivial solution to the gap equation, the normal 
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state is unstable. We wish now to determine whether 
our ansatz of Eq. ([5]) is stable. As discussed above, if 
there were a region with unstable modes (i.e. positive 
imaginary parts of poles), then this would lead to the 
existence of a true pole at real omega, at the boundary 
of the unstable region. Making use of the symmetries in 
Eq. (PT]) . for the condensed case, poles of the retarded 
Green's function correspond to solutions of 

Af (^, p)Af (-C., p)* - Rf{uj)K^i~Lur = 0. (32) 

Unfortunately this expression is not simple, and nu- 
merical evaluation would be necessary to trace the be- 
haviour of all zeros as a function of momentum. How- 
ever, in order to understand the stability, we can instead 
consider separately zeros of the real and imaginary parts 
of Eq. ((32|) . If zeros of these two parts coincide for some 
Pc , there is a real pole, and thus instability for |p| < pc- 
It is clear the imaginary part should have a zero at a; = 
(measuring frequency from the common oscillation fre- 
quency fis), as the imaginary part of Eq. ([5^ is an odd 
function of uj. This zero physically corresponds to the 
divergence of the distribution function at w = 0. Numer- 
ical investigation suggests that this is the only zero of the 
imaginary part. Thus, we are interested in zeros of the 
real part, evaluated at a; = 0, but arbitrary p. 

It is clear there is a zero at w = 0, p = 0, corresponding 
to the symmetry under global phase rotations, but being 
at p = this does not lead to instability. From this pole, 
or alternatively working directly from the definitions of 
in Eq. and the gap equation (fTi|) . one can show 

that K^{uj = 0,p = 0) = K^{uj = 0). Thus, writing 
A = Ke [kC{uj = 0, p = 0)] ; instability occurs if there is 
a non-zero p solution of: 



A 



1 p2 



-A' = 



1 p2 



-AA] = 



2 2mph / 4 27Tiph \ 2mph 

which will exist if and only if A > 0. 

Physically, this says that the Goldstone mode will be 
unstable for < |p| < Pc if the "static compressibility". 
Re [ATf (0,0)] > 0. In equihbrium, the expression for the 
component K^{0,0) is real and negative, but including 
pumping and decay, there are regions where solutions of 
the gap equation, Eq. ((T4)) exist but which are unstable. 
Since Re [if^(0, 0)] is the real part of the second deriva- 
tive of the action w.r.t ip{u} = 0,k = 0), it can also be 
seen as a derivative of the gap equation, thus unstable 
solutions are characterised by a non-linear susceptibility 
that increases as coherent field increases. 

As a result, there are ranges of the parameters k, 7, /i^ 
for which neither the normal state, nor the ansatz of 
Eq. ([8]) arc stable. We have not investigated what alter- 
nate stable solutions might exist under these conditions, 
however the existence of a real pole in the response at 
a non-zero momentum might suggest one should inves- 
tigate the possibility of a coherent field at non-zero p. 
Such a possibility would not be too surprising, as spon- 
taneous pattern formation is seen in laser systems with 
a continuum of modes^^. 



VI. NUMERICAL ANALYSIS OF THE 
MEAN-FIELD 

A. Phase diagram 

Having discussed the conditions under which the uni- 
form, single-frequency condensed solution is stable, we 
may now consider an effective phase boundary — i.e. 
find the ranges of parameters for which there is a stable 
condensed solution. For numerical analysis we choose all 
baths to be individually in thermal equilibrium. How- 
ever, as the baths need not be in equilibrium with each 
other, the system can still be far from thermal equilib- 
rium. Since the cavity photon modes start at energies 
much above the zero for bulk photon modes, we take the 
chemical potential of the decay bath to be large and neg- 
ative. In addition, since at room temperature the popu- 
lation of the bulk photon modes at the energy of cavity 
modes is negligible, we consider the decay bath to be al- 
ways at zero temperature. In the following we will first 
present calculations at zero pumping bath temperature, 
with a delta function density of states i.e. ea = e, and 
at zero detuning. Following that we will then analyse 
the influence of finite temperature of the pumping baths, 
and of inhomogcncous broadening of cxcitons. At zero 
bath temperature, the bath distributions arc entirely de- 
fined by their chemical potentials, and so there remain 
three control parameters, fj-B, 1, i^- Note that in this case 
the pump and decay baths are at the same temperature, 
but have very different chemical potentials, thus leading 
to a particle flux through the system, driving it out of 
equilibrium. In Fig. [51 we illustrate the boundary as a 
function of fis, 7, ^ by plotting its section in two planes; 
the plane of fixed k [Fig.[6l[a)], and the plane of fixed 
[Fig.Hb)]. 

It is worth noting that, for /is < 0, and fixed k,pb 
there is both an upper and lower critical 7. The maxi- 
mum 7 is always present (i.e. even if /.ts > 0), and results 
because increased coupling to the bath causes dephasing. 
Let us discuss the origin of the minimum critical 7. If the 
bath is at zero temperature, it pumps only that part of 
the effective excitonic density of states with energy less 
than the bath chemical potential /i^. If there is no in- 
homogeneous broadening (i.e. ea = e) then the effective 
exciton density of states is set entirely by its coupling to 
the baths; i.e. it is Lorentzian with width 7. Thus, the 
efficiency of pumping depends on how, by broadening the 
excitonic energy, the pump leads to a non-zero density of 
states below the chemical potential fiB- As a result, at 
7 = there is no pumping, and so no condensation, and 
a minimum 7 is required before there is sufficient gain 
to overcome the decay. If there is inhomogeneous broad- 
ening of exciton energies, or the pumping baths are at 
finite temperature, this effect is less significant, as is seen 
in Fig. [i 

From the boundaries of the stable region, it appears 
that a uniform condensed stable solution is only possible 
if K < kq, with kq ~ 0.2(7. The origin of this upper 
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FIG. 6: (Color online) Phase boundaries at zero tempera- 
ture, no inhomogeneous broadening. Panel (a): Fixed decay 
rate, k. Panel (b): Fixed chemical potential, \ib (Note that 
/is > implies inversion in the pumping bath but not nec- 
essarily in the system). The insets show in detail the region 
of small K, and 7. Solid lines mark the limit of stability of 
the normal state. Dashed lines mark the limit of stability be- 
tween the uniform condensed state, and some other unknown 
state. The uniform condensed state is stable in the shaded 
regions. The asterisk in panel (b) marks a point of fixed k, 7 
for comparision between the plots. 



critical k requires further investigation. 

B. Coherent fields and densities 

As well as the phase boundary, one may study the evo- 
lution of a number of properties of the condensate — e.g. 



trate both the regime of weak coupling to baths, where 
the results are similar to those in thermal equilibrium, 
and also strong decay and pumping, for which the results 
are instead comparable to the laser. For comparison, the 
value of ^ and the fermion-pair (excitonic) density in 
the normal state are shown, which connect smoothly to 
the condensed quantities, as expected for a second order 
phase transition. Note that for k = 0.15,7 = 0.9, 1.0 the 
excitonic density Pcxciton > 0.5 indicating inversion as is 
expected in the lasing case. 
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FIG. 7: (Color online) Properties of the system as a function 
of bath chemical potential. Plotted for T = 0, and = e. 
Panel (a): k = 0.02gi, panel (b) k = O.lSg, values of 7 as indi- 
cated. Within each panel, four graphs are shown. Top: The 
common oscillation frequency in the ansatz of Eq. ([HJ, mea- 
sured from a;p=o = 0; Second: Density of condensed photons; 
Third: exciton density from Eq. (|16fl : Bottom: Polarisation, 



mean-field density of condensed photons, excitonic gi^en by ^S^T^V'/- Solid lines indicate where a stable con 



density [from Eq. , and thus the total mean-field den- 
sity, being the sum of condensed photon and exciton den- 
sities, polarisation {a^h) (where | {a^b) p gives the num- 
ber of condensed fermion pairs - excitons), and common 
oscillation frequency jis- These are shown in Fig. [3 for 
two values of k and a range of different 7, chosen to illus- 



densed solution exists; dotted lines are the (unphysical) result 
of the uniform condensed solution when it is unstable. Dashed 
lines for ^,3 and pcxciton show the comparable quantities in the 
normal state. 
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C. Influence of bath's temperatures and excitonic 
density of states 

We now consider the effects of finite bath temperature, 
and of the inhomogeneous broadening of the exciton en- 
ergies. As such calculations are numerically intensive, 
we present a limited, but illustrative set of results. In 
Fig. [51 the equivalent of Fig. [D^b) is shown, but with 
a Gaussian density of states, and at small but non-zero 
temperature of the pumping bath (the decay bath, of 
bulk photon modes, is still at T=:0). One can clearly see 
that by adding inhomogeneous broadening, a,, = O.lSg, 
the lower critical 7 has been modified, and for large 
entirely eliminated. The inset of Fig. [5] shows a higher 
temperature, for which none of the curves show any lower 
critical 7. 
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bution is too. This means we can study the influence 
of dephasing due to pumping and decay separately from 
the influence of non-thermal distribution functions. This 
also allows direct comparison to the equilibrium limit, 
which, as discussed in Sec. IIV Al should be recovered as 
K — > 0, 7 0. This is illustrated in Fig. [SI where the 
critical bath temperature as a function of system density 
is plotted (and for comparison, the critical /is at each 
temperature is also shown). It is apparent that the pres- 
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FIG. 9: (Color online) Critical density (and associated crit- 
ical bath chemical potential) at a given non-zero bath tem- 
perature. Evaluated for a Gaussian density of states, with 
cte = 0.15g and values of n and 7 as indicated in the legend. 
The dotted line marks the limit k ^ 0, 7 ^ 0, for which 
the equilibrium result, with distributions set by the pumping 
bath is recovered. 



FIG. 8: (Color online) Phase boundary for constant chemical 
potential, /is, as in Fig. |6jb), but with a Gaussian distribu- 
tion of excitonic energies, = 0.15g, and non-zero tempera- 
ture (top row T = O.Olp, bottom row T = O.bg) 

. As a result, the requirement for a minimum coupling 
strength, 7, before a transition occurs is removed for some 
phase boundaries. Solid lines, dashed lines, and shaded 
region mark instability of normal state, instability of uniform 
condensed state, and stable condensed region as in Fig. [B] 

One can also plot a phase boundary at fixed 7, k as 
a function of pumping bath temperature T and /^s, or 
alternatively derive the excitonic density /Oexciton from 
Eq. to plot the boundary as a function of T and 
density. By doing this we can investigate the influence of 
decoherence and particle flux introduced by pumping and 
decay on the phase diagram, which can still be significant, 
even if the system distribution function would be close 
to thermal. For the parameters chosen for the figures, we 
are in the regime of densities where the phase transition 
is well described by mean-field theory, and so the number 
of incoherent photons at the transition is small. Thus, in 
this regime, the distribution function of excitons below 
and at the transition is set by the pumping bath; thus 
if the pumping bath is thermal, then the exciton distri- 



ence of pumping and decay shifts the phase boundary 
to higher densities, and that the p as T be- 
haviour seen in equilibrium does not survive. Physically, 
this increase of critical density is due to the decoherence 
introduced by pumping and decay. The behaviour at 
T ^ is unsurprising, as the limit p corresponds 
to the equilibrium chemical potential fi ~^ —00. In the 
presence of non-zero decay rate k, one requires a non- 
zero effective gain (imaginary part of gap equation) , and 
so no solution exists with fis —00 even at T = 0, i.e 
the critical density never goes to zero. 



VII. FLUCTUATIONS IN CONDENSED STATE 
TO ALL ORDERS IN PHASE 

The low energy modes of the broken symmetry sys- 
tem correspond to slow phase variations. Since there 
is no cost to global phase rotations, the action depends 
only on derivatives of the phase, and so phase fluctua- 
tions may become large. Thus, describing '4> = ipQ + Sip 
and considering only terms to second order in Sip may 
underestimate how phase fluctuations reduce long range 
coherence. Therefore, we will instead consider the pa- 
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ramcterisation -0 = ^ pq + Tre*"^, and evaluate correlation 
functions of i}} in terms of the correlation functions of 
amplitude tt and phase 0, including the phase fluctua- 
tions to all orders. In equilibrium, the effect of phase 
fluctuations on the field-field correlator is responsible for 
the reduction from long range order to power law corre- 
lations in two dimensions, and so has been much studied 
(see e.g. Refs. HIsQl . Here, in order to calculate the 
luminescence and absorption spectrum, we will however 
need also to include density fluctuations. 

Combining such a reparameterisation of the fields 
with the non-equilibrium Keldysh formalism requires a 
little care. The first important consideration is that 
the parameterisation requires one to work with fields 
where is macroscopic. This means wc should re- 

parameterise the fields i/'/ j V'fe defined on the forward and 
backward contour (see Sec. IIIip . as opposed to the fields 
ipq,ipci, since (l^gP) is not macroscopic. This consid- 
eration is similar to the fact that the parameterisation 
should be done for the fields as functions of space and 
time rather than functions of p and uj. The second con- 
sideration is that, in calculating the physical correlation 
functions, 2?^'^, this will involve cross terms between 
the two branches, and so one must keep track of which 
branch tt and (j) are on. 

The technical details of how to derive the field-field 
correlation functions in terms of amplitude and phase 
Green's functions are presented in Appendix |B] For 
the forward Green's function (corresponding to lumines- 
cence), the result is found to be: 



1 



^ Uv<^{Q, 0) + tV<^{0, 0) - ^V<Jt, r) - iV<^{t, r)] ' 



{- [*25|0(O,O)-zP<^(i,r)]} (33) 



exp • 



The above procedure includes amplitude fluctuations tt 
and gradients of phase fluctuations V0, dt(j} to second or- 
der as they both have restoring force, and cost energy, so 
that they are expected to be small. The phase fluctua- 
tions 4> however may be large and in the above result are 
taken to all orders. 

To calculate the luminescence spectrum one must then 
Fourier transform the result V^^^lt^r) to give the spec- 
trum in frequency and momentum space. The first term 
in the braces in Eq. (|33[) proportional to po describes the 
emission from the condensate which is now broadened by 
the exponential term containing the phase fluctuations. 
It is clear that the phase fluctuations determine the con- 
densate lineshape and the decay of spatial and temporal 
coherence. An example of luminescence as given by Eq. 
([55]) is shown in Fig. (TUl We will discuss its features in 
Section EIEl 
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FIG. 10: (Color online) Photoluminescence from the region 
of small LJ and small p (shown it terms of an angle of emis- 
sion 6 = tan~^(cp/ci;o), calculated including effects of phase 
fluctuations to all orders. Both panels have n = 0.02(7, 7 = 
0.2(7, hb = O.Og. The upper panel has T = O.lg as in the mid- 
dle row of Fig. 111! while the lower has T = O.Olg, for which 
the features described in the text appear more sharply. 



If one were to assume phase fluctuations were small, 
then this expression could be expanded to linear order in 
Green's functions, and one would find: 



^T^^^i,(t^r)=po{l 



*2?^.(0,0) 



*2?1(0,0) 



4/90 



This is instructive, as the second line describes the fluc- 
tuation Green's function iT>f^fg^{t,r), obtained taking 
the fluctuation fields to second order, while the first cor- 
responds to a depleted condensate density. Such a lin- 
earisation would describe the luminescence as a sum of 
two terms; a condensate term, which due to its lack of 
space or time dependence would be a sharp peak, and 
a fluctuation term. Furthermore, if one were to consider 
the frequency spectrum of fluctuations by integrating this 
linearised form over momentum one would have a sim- 
ple power law form, with a power depending only on 
the dimension^^, and not on parameters of the system. 
By allowing phase fluctuations to be large, and keeping 
the phase-phase Green's function in the exponent, the 
condensate acquires a lineshape as a result of phase flue- 
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FIG. 11: (Color online) Spectral weight, photoluminescence and absorption spectra, as a function of emission angle, 
tan~^(cp/ci;o)- For all graphs, k = 0.02g and T — O.lg. Top row: Uncondensed case, 7 = 0.2g,^_B ~ —0.5g. (cf parame- 
ters in Fig.Oand Fig. 2} Middle row: Condensed case, 7 — 0.2g,iiB = O.Og. Bottom row: Condensed case, 7 — 0.5g,fj.B = O.Og 
(transition to weak coupling). 



tuations, and this lincshape can in the equihbrium hmit 
recover the standard power law correlations seen in two 
dimensions. The form of this lineshape is discussed fur- 
ther in Sec. IVIIAI 

However, for uj,p far from zero, such linearisation does 
not introduce any major changes; the effects of large 
phase fluctuations matter mostly at large times. Large 
fluctuations between fields separated by small t or r 
would imply large gradients, and thus have a large energy 
cost. Thus, Fig. [11] illustrates the absorption, lumines- 
cence and spectral weight over large ranges oi uj,p using 
a linearised approach [which at this large scale coincides 
with the full expression given by Eq. ([551) ] while Fig. [TUl 
obtained from the full expression of Eq. ([55]) . shows the 
effect of phase fluctuations at small 

For the detailed analysis of the features of the lumi- 
nescence spectra we refer to Ref. [H. Note that for large 
Lu,p as shown in Fig. [TT] the main features of the non- 
equilibrium spectra are similar to those predicted for 
equilibrium condensation in Refs. [30ll43l . In the nor- 



mal state one can see the upper and lower polariton 
modes (top row of Fig. [TT|) in the spectral weight and 
absorption, and only the lower polariton in the lumines- 
cence as the upper polariton is not occupied at this low 
power. When system condenses (middle row of Fig. [TT|) 
the structure of modes changes dramatically showing the 
pairs of phase and amplitude modes above and below 
the chemical potential. Finally, when the coupling to 
the pump baths; i.e. the pumping strength is further in- 
creased (bottom row of Fig. [TT|) the system crosses to 
weak-coupling regime and the polariton splitting is sup- 
pressed. In Fig. [11] the occupation of the excited states 
will not be thermal, in contrast to the analogous figures 
m Refs. MM 

this is however not easy to observe on 
these contour plots. Also since Fig. [11] corresponds to 
pumping baths at finite temperature in contrast to zero 
temperature in Ref. [l^ the sharp occupation edge visible 
there is here smeared out. However the main qualitative 
difference between the spectra of a pumped decaying con- 
densate presented here and that for a closed system given 
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Refs. MM 

is most visible on small u;,p scale as pre- 
sented in Fig. [ini This will be discussed in detail in the 
Section IWAl 



and thus depend on r and ct. For small p, the lead- 
ing term in the expansion for both cos(cpt) and Joijpr) is 
quadratic, and so the lower cutoff for the integral is given 

by p ~ l/^/r^ + cH^. Thus, 



A. Condensate Lineshape - effects of dissipation 
and low-dimensionality on decay of correlations 

The long range field-field correlations are influenced by 
the properties of the soft phase modes; i.e. the Goldstone 
or Bogoliubov mode o^^i^^ . By considering the asymp- 
totic behaviour of the phase-phase correlator at small fre- 
quencies and momenta, one can thus find the asymptotic 
form of the field-field correlator. In an equilibrium two- 
dimensional system, the long distance field-field correla- 
tions decay with a power law below the BKT transition. 
We will now investigate how this asymptotic behaviour 
is affected by the presence of pump and decay. For con- 
venience let us rewrite Eq. (|B4p . assuming an isotropic 
system: 

iV<,^{t, r) = po [1 + O (1/po)] exp [~f{t, r)] , (34) 

(35) 

Here Joipr) is a Bessel function, from the integration over 
azimuthal angle. We are thus interested in the limits 
f{t = 0, r oo) and f{t oo,r = 0), describing the 
large distance and long time decay. 

For comparison, let us first summarise how this method 
reproduces the standard result in the equilibrium case. 
In equilibrium, the distribution function is a constant 
matrix F{uj) = 2nB{io) + 1, and so: 

iV<^{Lu,p) = 1 (F(^) - 1) {iV^^{u,p) - iV^^{u,p)) 



= jiB(c^)(-2)Im [V^^iuj.p)] 



(36) 



For an equilibrium coherent system, the low energy 
modes will be the linear Goldstone modes of the form 
Lo = cp. By analytic continuation of the imaginary time 
(Matsubara) Green's function, one finds: 



Im[l?« (w,p)] =Im 



-C 



= -^{S{u;-cp)-6{u; + cp)) (37) 

And so, combining Eq. ^5^, Eq. and Eq. one 
finds that the singular contribution to f{t, r) is given by: 



fit,r) 



C 



1//3 



"dp 

2tt(3c Jo p 



1 — Jo{pr) cos{cpt)] 



(38) 

where /3 is inverse temperature, and the effect of the 
thermal distribution has been approximated by the up- 
per cutoff of the integral. The lower cutoff is con- 
trolled by how Jo(pr) cos(cpt) approaches 1 as p ^ 0, 



f{t,r) ~ ?7ln 



]3c 



Thus, one recovers the standard result, and logarithmic 
behaviour of f(t, r) leads to power decay of correlation 
functions, with i] oc ksT/pQ. One can further use this 
result to find the form of the peak in the luminescence 
spectrum, C{uj,p) cx (c^p^ -|- w^)^''"^^''^, and the inte- 
grated luminescence (i.e. angular profile^) N{p) cxp''"^. 

Let us now consider the asymptotic form of the Green's 
function in the non-equilibrium case. We shall first con- 
sider the retarded Green's function, as the poles of this 
function describe the normal modes; the result of calcu- 
lating will, as discussed later, be to introduce the 
population of these modes. The retarded Green's func- 
tion, using the notation of Eq. (|20p can be written as: 



C 



As discussed in Sec. IV CI the gap equation implies that 
Ki{lj = 0,p = 0) = A'^(ti; ~ 0). Combining this with 
the symmetries in Eq. ((2T|) . one can show that the most 
general expression, to quadratic order in p, lo in the de- 
nominator can be written as: 



C 



c?p'^ + 2iujx 



(39) 



where C, c and x are coefficients to be derived from the 
full expressions. Without pumping and decay, a; = 0"*", 
and one recovers the equilibrium result. With non-zero 
X, the poles of the Green's function, which define the low 
energy modes of the system, have the form 



c^p'^, 



and are thus diffusive, rather than dispersive for p < 
x/o^. This can be clearly seen in the luminescence shown 
in Fig llOl At low momentum, where the real part of the 
pole vanishes, but the imaginary part does not, the lu- 
minescence is dispersionless (i.e. flat), but broadened. 
Such a form should be generic for broken symmetry in 
a pumped decaying system, and indeed the same form 
has been recently seen in a related context, of coherently 
pumped polaritons in photonic wires, described as an op- 
tical parametric oscillator—, as well as in a more generic 
model^. This result also shows why it was so important 
to have solved a complex gap equation, rather than just 
adding decay rates to the equilibrium model. Adding 
phenomenological decay rates "by hand" would lead to a 
form of the retarded Green's function: 



C 
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Such a form does not describe a system with sponta- 
neously broken symmetry, as there is no pole at w = 
0,p = 0, and thus such an approach misses the appear- 
ance of a diffusive mode. 

Let us now consider 2?^^, and thus the effect of the dis- 
tribution function. As was discussed in Sec. lV A[ the dis- 
tribution function can be expected to diverge at the en- 
ergy where the imaginary part of the denominator of the 
retarded Green's function vanishes. This is clear atuj~Q 
(measured relative to the common oscillation frequency 
fis), due to the presence of a real pole at a; = 0,p = 0. 
However, this divergence will be exactly canceled by the 
vanishing of 'D^{uj,p) — 'D^{uj,p) as ui — *■ 0, since both 
the divergence and the vanishing are due to the same 
imaginary part. Thus, near w = 0, the asymptotic form 
of Pf.(w,p) is the same as that of \'D?-Juj,p)\'^, i.e.: 



C2 



(cj2 — c?p'^Y + 4a>2a;2 



The effect of the distribution will be to introduce some 
upper energy cutoff. Thus, the equivalent of Eq. ([35]) is: 

/(^'0 = T^/ -[l'Mpr)d{p,t)] (40) 
zc X Jq p 

where the time dependence is described by: 



d{p, t) = e" 



\J — (?p^ 



sinh I ya;^ — (p-p^t 



+ cosh ( \/ — (P'p'^t 



(41) 



Equation (|40p has a similar interpretation to Eq. 
a large p cutoff from the distribution function and a short 
distance cutoff set by the coordinates. For a thermal 
distribution function, the upper cutoff would be given 
by 1/^c — fcsT/c. Although the photon distribution 
in the pumped decaying system is not thermal, if the 
pumping and decay baths are thermal (as considered ear- 
lier), then the photon distribution will vanish for large 
enough energies. As such, we will write 1/^c — ^'max/c, 
where iSmax depends on both pumping and decay, and 
would reduce to fcsT in equilibrium. The result is thus 
/(i, r) ~ ry' ln(l/Q^c)j where Q is the lower cutoff. How- 
ever, the form of the lower cutoff can be different, and 
depends on the relative values of r, ct and cjx. In the 
two regions of interest defined at the start of this section, 
one finds: 



— 1= if 7' ~ 0, < ^ cx), 
i if r ^ OO, t ~ 0. 



Inserting this cutoff, one finds 



[t]' l2)\ii{cHlxil) ifr: 
77' ln(r/^c) if T 



0, t ^ 
» 00, i ■ 



00, 
- 0. 



(42) 



(43) 



Thus, there is still power law decay, but due to pumping 
and decay the powers for temporal and spatial decay do 
not match, and since 77' may depend on x, both power 
laws will differ from equilibrium. 

Since the long time decay is power law, the lineshape 
will also have a power law divergence at low frequency, 
and as such there is no well defined condensate linewidth 
in an infinite system. In fewer than two dimensions, i.e. 
in a ID system^!, or a fully confined system such as a 
laser with discrete modes, the long time decay will be 
exponential, and so a linewidth can be found in such 
systems. The crossover between power law and exponen- 
tial decay in a large but finite 2D system is discussed 
in Sec. IVIIll In three dimensions, the limit of f{t,r) at 
large times and distances is finite (as opposed to diver- 
gent as in two, one or zero dimensions). As a result, there 
is phase coherence to arbitrarily large distances, and so, 
writing the asymptotic values of f{t,r) as /oo there is a 
contribution to the luminescence that goes like: 

iV<,^{u;,p) = JdtJ dhpoe-f-e^'^'+^P-'^ + ... 

i.e., in an infinite homogeneous 3D system, there would 
be a peak at w 0, p = 0, with a peak height given 
by the condensate density, which is depleted by phase 
fluctuations. 



VIII. FINITE SIZE EFFECTS 

In the previous section we discussed how the contin- 
uum of phase modes leads, in two dimensions, to loga- 
rithmic phase-phase correlation functions as a function 
of distance and time. In this section, we consider how 
confinement, which leads to a discrete spectrum of phase 
modes will modify that result. In a confined system, 
there will not be translational invariance, and so the field- 
field correlation function will in general depend on both 
positions, rather than just on separation. However, if 
we are interested in the equal-position, long-time limit, 
which is relevant for the lineshape, we can then write: 



dwC\Mr)\H^ 



27r (tj2_ (2)2 4, 4^23.2 



where we have introduced the wavefunction y:>n{r) and 
energy („ of the n"^ phase mode. It is clear that if 
V'ni'i") = e'P"'", and C„ = cp„, we recover the previous 
result. 

Let us now discuss briefly the energy spacing A of 
phase modes Schematically, for a box of size R, one 
has A = c/R, i.e. the sound modes, with discrete mo- 
mentum spacing. In contrast, the energy spacing of sin- 
gle particle states in such a box would he S = \/2mR^. 
Since the sound velocity increases as condensate density 
increases, one can have A 3> (5. [NB in a harmonic trap. 
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the Thomas-Fermi radius and the sound velocity have the 
same dependence on pq, so the phase mode level spacing 
is the single particle spacing^. A harmonic trap is how- 
ever a special case in this regard.] 

To understand how discrete mode spacing modifies 
f{t,r,r), let us first reconsider how the logarithm term 
arose from the integral. Schematically, we had: 



dp 



^Z^' dp 
P 



— = ^ + In 



1 



i.e., the dependence on the coordinates, via the cutoff Q is 
logarithmic, as the contribution from p < Q is constant. 
For the discrete sum, after integrating over u, instead of 
Eq. (gOl) we have: 



fit,r,r) 



'2x 



N 

E 



(44) 



To summarise, if temperature is low enough (or in the 
case of non-thermal distribution the relevant energy to 
which the modes are occupied is small enough), phase 
fluctuations are frozen out, as one expects. If phase fluc- 
tuations are not frozen out, there are two limits; at long 
enough times, one always sees linear growth of fluctua- 
tions, resulting in the exponential decay of field-field cor- 
relations, and recovery of the standard laser lineshape^. 
However, for large enough systems, so level spacing is 
small, there is a range of time delays during which the 
growth of phase fluctuations is logarithmic in time, giv- 
ing rise to a power-law decay of fleld-fleld correlations, as 
one would expect in the infinite system. 



A. Self-phase modulation 



with d{p,t) as in Eq. ([^T|) . The upper cutoff is in- 
troduced here by truncating the sum at N such that 
Cn = c/^c — £'max- Considering the long time limit, 
this sum can also be split into two parts; for modes 
(^n <C ly^x/t the summand is effective energy indepen- 
dent, while for C„ ^ \/x/t, with the density of states in 
2D, one recovers a log divergence. However, the existence 
of these two parts depends on t he relative values of the 
energy of the lower cutoff \/x/t, the upper cutoff i^max, 
and the level spacing A. We assume E'max ^ \/x/t, 
which just means considering long enough time delays, 
and so there are three important cases: 

1. A <C \/x/t <C i^max- In this case there are many 
terms contributing to both the small and large Cn 
sums, and so the the result is as for the integral: 
schematically f{t,r,r) = 1 + ln(£'„iax\/t/a;), and 
there are power laws, as in the infinite system. 
This case cannot however persist to arbitrarily large 
times. 

2. \pxjt < A ^ £'max. At long enough times, the 
previous case will switch to this case. Here, there 
are only a few terms in the low energy contribution. 
A characteristic term, for C,n x gives d{C,n/c, t) ~ 
1 — C,^t/2x. Since the number of low energy modes 
is now of order 1, rather than of order x/tA'^, the 
contribution from these modes is of order t/x, and 
not of order 1. Thus, the dominant contribution 
is f{t,r,r) ~ {-kC /2x){t/2x), and so the decay of 
field-field correlations is exponential as in a single 
mode case. 

3. yjx/t <^ Eraax ^ A. In this case, no phase 
fluctuations are populated, i.e. no terms survive 
in the sum, and so the entire system is coher- 
ent. Using A = c/i?, this condition is equivalently 
i? <C = c/i?niax, i-c. the "thermal length" is 
larger than the system size^. 



The analysis so far shows how, due to finite size, the 
power law correlations associated with a continuum of 
modes change to the exponential decay of correlations as- 
sociated with phase diffusion of a single mode. There has 
been previous work on extending the picture of phase dif- 
fusion of a single mode due to pumping noise*^ to the case 
of interacting systems, for which there is an additional 
source of noise from self-phase modulation (SPM)^i2ii^. 
These works suggest that the phase decay rate can be 
written as a; ~ (Fq -|- p^Xspm)/ Po, where Fq is the noise 
due to pumping, po the condensate density, and Xspm 
proportional to interaction strength. We wish here to 
comment briefly on the origin of the SPM term, and how 
it may be modified in the case of many interacting modes, 
with respect to the case of phase diffusion of a single 
mode. 

The presence of a SPM term can be understood by 
considering the evolution of a coherent state, e^^"^'^ |0) = 

J2niV p"/"OI"-)- Foi' ^n interacting single mode system, 
the number states are eigenstates, and have energies like 
En = an + bn^, thus different number states evolve at 
different frequencies, and mutually dcphasc, leading to a 
dephasing rate Xspm — bpo- Thus, SPM occurs because 
number states, not coherent states are eigenstates of the 
single mode Hamiltonian. The eigenstates of the many 
mode system, including coherent interactions between 
modes, such as ^\^\'ipp^_p are neither number states 
nor coherent states, but are instead better described by 
Noizerres-Bogoliubov states^. Such states are superpo- 
sitions of terms with different divisions of particles be- 
tween the condensate and non-condensed modes; while 
they may be eigenstates of total number, they are not 
eigenstates of the number of particles in a given mode, 
and they lower energy because of the coherence between 
the different modes^. As such, when considering systems 
with a continuum of interacting modes, it is not clear that 
SPM terms should exist, or if they exist, should have the 
same form. 
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IX. CONCLUSIONS 

In conclusion, wc have studied steady-state sponta- 
neous quantum condensation in a non-equilibrium Bose- 
Fermi system with pumping and decay, and consequent 
flux of particles. In order to study the effect of large 
phase fluctuations in the broken symmetry system, it was 
necessary to extend the path-integral Keldysh formal- 
ism to deal with a reparameterisation in terms of phase 
and amplitude fluctuations, for fields on the forward and 
backward time contours. Wc have shown that the mean- 
field properties of a pumped and decaying condensate 
can be described by a complex analogue of the Gross- 
Pitaevskii equation in the BEG regime (or equivalently 
the gap equation in the BGS regime). The real part of 
this self-consistency equation relates the coherent field 
to the system's non-linear susceptibility as in the case of 
equilibrium condensation, while the imaginary part re- 
flects how the gain and decay are balanced, as in a laser. 
We further show that it is crucial to satisfy this complex 
self-consistent equation in order to get the correct collec- 
tive mode structure, reflecting the broken symmetry. 

We have analysed the solutions of this complex gap 
equation and examined their stability. Surprisingly, de- 
spite non-thermal distributions, the instability of the nor- 
mal state is analogous to that in thermal equilibrium, 
where the normal state becomes unstable when the chem- 
ical potential, at which the Bose-Einstein distribution di- 
verges, reaches the bottom of the system's spectrum. In 
the non-equilibrium case, the system's distribution, al- 
though far from thermal, develops a divergence at some 
energy. When, by tuning parameters of the system, this 
energy is brought to coincide with an effective pole of 
the system's Green's function, then the normal state be- 
comes dynamically unstable and the condensation tran- 
sition takes place. We have also shown that whenever 
there is a condensed solution, the normal state becomes 
dynamically unstable, and so there is no ambiguity as to 
which state the system would choose. However, we have 
found a range of parameters where both the normal and 
the uniform harmonic condensed solutions are unstable, 
suggesting either more exotic, perhaps chaotic, dynamics 
or spatial pattern formation. 

We have analysed the non-equilibrium phase diagram 
as a function of the decay and pump parameters, and 
have found both the low density condensed solutions 
when pump and decay strengths are relatively small, 
as well as the high density, inverted, laser-like solutions 
when pump and decay are comparable to the inter par- 
ticle interactions. When applied to microcavity polari- 
tons these regimes reflect the spontaneous condensation 
of strongly coupled photon-exciton modes at relatively 
small pump and decay powers, and the crossover to the 
weak-coupling regime and the photon laser at large pump 
powers. It is important to stress that even if the system 
distribution is close to thermal the presence of pump and 
decay, i.e. particle flux, result in a higher critical density 
at a given temperature than in a closed system, and there 



is a non-zero minimum critical density even at zero tem- 
perature. 

Having analysed the fluctuation spectra and collective 
modes, we have found an important difference between 
condensation in a dissipative environment and that in 
closed systems: Although there is a real pole (undamped 
mode) at zero frequency and momentum, indicating bro- 
ken symmetry, the usual linear dispersion of the sound 
mode (Bogoliubov, Goldstone mode) at small momenta is 
now replaced by diffusive behaviour (i.e. a broadened but 
flat dispersion) ; this questions the possibility of superflu- 
idity on large time and distance scales. This qualitatively 
new structure of the collective modes is visible in the 
luminescence and absorption spectra, and it affects the 
field-field correlations, i.e. decay of spatial and temporal 
coherence, and the condensate lineshape. For example, 
in the 2D system dissipation changes the usual power-law 
decay of spatial and temporal coherence, replacing it by 
one where the powers for temporal and spatial decay do 
not match. 

It is instructive to place our treatment of non- 
equilibrium quantum condensation in the context of 
other works on dynamic effects in polariton systems. 
Much of the literature concentrates on Boltzmann-like 
rate equations. Such an approach allows one to study 
the effect of pumping and decay on the occupation of 
modeSf^iSi^ but is not able to account for the changes 
to the excitation spectrum and the density of states 
(which are particularly dramatic as the system crosses 
the phase transition) due to the pumping, decay and 
presence of the coherent field. In contrast, field the- 
oretical studies presented here selfconsistently account 
both for arbitrarily large changes to the excitation spec- 
trum as well as changes to the occupation of this spec- 
trum. Such approaches are thus well placed to study 
the phase transition between the non-condensed and con- 
densed states, and in addition the crossover between 
strong and weak coupling regimes. A closer approach 
to the field-theoretical approach presented here would be 
the evolution of the off diagonal parts of the full density 
matri»22i^. It was however only recently that qualitative 
changes to the spectrum have been calculated using the 
density matrix approach (in the context of parametric 
emission from photonic wires) in Refs. [2^l25l . A further 
distinction is between single-mode models, in which one 
expects phase diffusion (e.g. Refs. [S^IH) and exponen- 
tial decay of correlations as in lasers, and models with 
a continuum of modes, such as Refs. [T8ll2^l25l and this 
paper. 

Finally, in this paper wc have analysed how the finite 
size of the system affects the decay of temporal coher- 
ence. This is particularly important for the understand- 
ing of recent experimentsi^, as well as for providing a 
connection to similar analysis for single mode photon 
lasers, which are still used as the basis to describe the 
decay of coherence in atom and polariton lasers. The 
key difference between the output from the condensate 
and from a single mode laser is that in the condensate 
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there is a continuum of modes, and so spatial fluctua- 
tions play an important role — in 2D they destroy the 
long-range order and lead to a power-law decay of corre- 
lations. Including such spatial fluctuations, the growth 
of phase fluctuations as a function of time is logarith- 
mic, which gives power-law decay of temporal coherence, 
rather than the exponential decay expected for a single 
mode. In single mode systems such as the laser, there 
are no spatial fluctuations, and so the decay of coher- 
ence is determined entirely from the phase diffusion of 
this single mode. However, if one takes a continuum sys- 
tem, and reduces its size, the energy spacing of modes 
becomes larger, and so the number of modes whose ener- 
gies are low enough to be relevant decreases, eventually 
recovering the single mode limit. We have identified two 
regimes in the finite system: Where the level spacing is 
larger than temperature, and so spatial fluctuations are 
essentially frozen out, resulting in an exponential decay of 
correlations as in a single mode laser; and where the level 
spacing is small with respect to temperature, so one gets 
a power-law decay of temporal coherence at short times 
as in the infinite system, crossing over to exponential de- 
cay at larger times. 

The qualitative implications of our results are general, 
and can apply to any BEC or BCS condensate which 
is subject to dissipation. The immediate applications 
of this analysis are for polariton BEC, which are natu- 
rally faced with significant pumping and decay processes. 
However, the techniques and results developed here, can 
be of use in understanding a wider class of broken symme- 
try dissipative systems; for example resonant parametric 
oscillators, and atom lasers, where coherence, dephasing, 
and the interaction of many modes are all relevant. 
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The expression for fermion-pair (exciton) density (fT5 
at r = reduces to: 
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APPENDIX A: GAP EQUATION AT T=0 

In the limit of T = 0, the integrals in the gap equation, 
Eq. (|14p can be evaluated in terms of elementary func- 
tions. This makes numerical analysis of the equations 
much easier in this limit. The results of this analysis are 
presented in Sec. IVIl for completeness we show the ex- 
plicit expressions at T = here. At T = 0, the bath 
distributions take a simple form Fb{uj) = sign([j — /Ib) 
and Fa{Lu) = sign(ti; + /jb) and so the real and the imag- 
inary parts of the gap equation become: 
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APPENDIX B: EVALUATION OF FIELD 
CORRELATIONS IN TERMS OF AMPLITUDE 
AND PHASE FLUCTUATIONS 



To illustrate the idea of using phase and amplitude 
fluctuations, we will first present the simpler case of 
2?"'"'^, for which both fields are on the same branch, and 
so we may drop all labels identifying which branch or 
Green's function (T or T) we are considering. Then, writ- 
ing ± for the time, and coordinate indices (T±i/2,R± 
r/2), one may write: 



The square root may be expanded to second order in the 
density fluctuations (as density fluctuations, unlike phase 
fluctuations, have a restoring force), thus: 
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iT^-Ji^^ — Pq 
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Introducing a current J. one may write the correlators in terms of a generating functional as: 
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By integrating over the photon field, the generating func- 
tional, Z[J] ~ (exp[. . . ]), can be expressed in terms of the 
correlators of amplitude and phase fluctuations. Defining 



J(^>P) - ( 2sin[(cjt + p-r) /2] 



one may then write 



Z[J] = exp 
where 
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4"!' 



Note that we use the standard definition of Green's func- 
tions so that iVab — (ab). 

To determine these correlators one may either recalcu- 
late the effective action by writing ip in terms of tt and 
4> and expanding to second order in tt and derivatives 
of (f>; or one may use the fact that at second order, the 



J 



amplitude-phase variables can be considered as a linear 
transform of S-ip and Sip, i.e.: 



[Sip 
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1 / 2po 2pa 

2Jp^ V « 



Note that this rotation relates the effective action ex- 
pressed in terms of these variables, and is not to be used 
in finding the final correlation functions I?^,t^- Thus, 
one can express the amplitude-phase Green's functions 
in terms of the Sip, Sip^ Green's functions as: 



f,R/A/K ^ L-D^/A/K ^ 



(B2) 



The T, T Green's functions can then be found from the 
retarded, advanced and Keldysh components by using 
Eq. ^ and: 



V 



,T,T 



1 



(B3) 



Thus, one may write the T or T field correlation function 
in terms of the phase and amplitude Green's functions: 



iT^i,^^ = Pa 



< 1 — > sm {ujt -(- p • r) — > [1 ~ cos [ujt -I- p • rjj 
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\\^^^^ cos(cjt-Fp • r)] '^-/-^(^'P) j I exp J - ^ [1 - cos (cjt p • v)]iV^^{uj,p) 



We can now address how to generalise this calcula- tion; the backward (absorption) will follow by swapping 
tion when the two fields are on different branches. We labels. Thus, repeating the above discussion, but keeping 
will consider the forward (luminescence) Green's func- subscripts on the fields one has: 
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1 + 



2po 



exp 



Then, as before, introducing a current, we may write this in terms of a generating functional. However, to keep track 
of labels, we shall need two currents, Jf and J;,, thus: 
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The calculation proceeds as before, but now the gen- 
erating functional is written in terms of: 
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J(w,_p) 
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There is thus an additional 2x2 structure of the Green's 
fimctions associated with branch labels, i.e. in block no- 



tation, 



With such an extended matrix structure, one can gen- 
eralise Eq. (|B1[) . and thus find the result: 



J 



E 

,uj,p 



This expression can be slightly simplified, since as ex- 
plained in Appendix [Cl below, one has that: 



= 0. 



(B5) 



LJ,p 



Using this result, and performing the Fourier transforms 
that appear in Eq. (|B4p . one then finds the final form for 
the Green's function, as given in Eq. 



APPENDIX C: ANALYTIC PROPERTIES OF 
GREEN'S FUNCTIONS 

Because the use of phase and amplitude variables 
forces one to work in terms of the physical Green's func- 



exp J - E (1 - exp [tiiut + p • r)])iV<^iuj,p) . (B4) 



UJ,p 



r 



tions, , iD^ , , and iV^ , it is necessary to con- 
sider the analytic properties of these Green's functions. 
As discussed in Ref. these Green's functions are not 
independent, but for t ^ one has 



(Cl) 



This lack of independence is implicit in Eq. 
and (|B3|) . repeated here for convenience: 
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1 
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(C2) 



However, at t = 0, Eq. (jCip does not hold. For the case 
of the field-field correlations, as discussed in Ref. the 
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correct regularisation leads to 



(C3) 



where TV is total particle number, and V is volume. The 
difference of form here is expected, as it encodes impor- 
tant information about the equal time commutation re- 
lations, 



The amplitude-phase Green's functions are found in 
the retarded, advanced and Keldysh basis, but to de- 
rive the field-field correlators we must rotate them to the 
forward and backward basis. Since this includes Green's 
functions of non-commuting operators evaluated at t = 0, 
such as Eq. (|B5p . it is important to reconcile Eq. ()C3|) 
with Eq. (jC2|) . Naively, such a reconciliation does not 
seem possible, however the resolution is that one must 
write: 



lim 



(C4) 

Thus, as one expects in a path integral formulation, op- 
erator ordering has been encoded via time ordering^. 
Written in terms of Green's functions as functions of fre- 
quency and momentum, the left hand side of Eq. (|C4[) 
would involve a conditionally convergent sum of terms 
that go like \/uj. Preservation of commutation relations 
thus requires correct regularisation of such conditionally 
convergent sums. The relations for the amplitude-phase 
correlation functions can be similarly found to corre- 
spond to the definition 



0) 



--{N + l)/V 



(t = 0) 



N/V 
(C5) 



With such a convention, the correct regularisation of 
the sum in Eq. (jB5P is then clear: 



iVl^it ^ 0+, r 0) - iVl,{t ^ 0-,r = 0) 



= 0. 



(C6) 
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